Low-energy effective theory and two distinct critical phases in a spin-^ frustrated 

three-leg spin tube 
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Motivated by the crystal structures of [(CuCi2tachH)3Cl]Ci2 and Ca3Co20e, we develop a low- 
energy effective theory using the bosonization technique for a spin- 1/2 frustrated three- leg spin 
tube with trigonal prism units in two limit cases. The features obtained with the effective theory 
are numerically elucidated by the density matrix renormalization group method. Three different 
quantum phases in the ground state of the system, say, one gapped dimerized phase and two distinct 
gapless phases, are identified, where the two gapless phases are found to have the conformal central 
charge c = 1 and 3/2, respectively. Spin gaps, spin and dimer correlation functions, and the 
entanglement entropy are obtained. In particular, it is disclosed that the critical phase with c = 3/2 
is the consequence of spin frustrations, which might belong to SU(2)fc=2 Wess-Zumino-Witten- 
Novikov universality class, and is induced by the twist term in the bosonized Hamiltonian density. 

PACS numbers: 75.10. Jm,75.40.Mg,75.50.Ee 



I. INTRODUCTION 



With the discovery of several spin tube materials, 
such as [(CuClatachrQaCyCla 1 , Na 2 V 3 7 2 and Cu 2 Cl 4 - 
D 8 C4S02 3 , etc., much effort has been made to investi- 
gate the low-lying excitations 4-9 , long range order 10 ' 11 , 
phase transitions 5,12 , magnetization plateau 13-15 and 
other ground state properties of spin tubes due to the 
inter-chain frustrated couplings. A spin tube, geometri- 
cally, can be seen as a multi-leg spin ladder with peri- 
odic condition along the rung direction (See Fig. 1 for 
the case of three- leg ladders). Therefore, the antifcrro- 
magnctic (AF) inter-chain couplings in the odd-leg spin 
tubes usually give rise to geometrical frustrations. For 
the typical spin tubes 4-6 [Fig. 1 (a)], the excitation is 
gapped when the inter-chain couplings arc identical, and 
becomes gapless when the interactions are different to 
some certain extent from each other. The opening and 
closing of the gap correspond to a dimerized phase and 
a critical phase with central charge c = 1, respectively. 
Recently, a few new spin tubes with more complex inter- 
chain couplings 7,12,15,16 [Fig. 1 (b) and (c))] have at- 
tracted much attention, in which frustrated inter-chain 
couplings could generate twist terms that may lead to 
unknown quantum phases in the ground state of the sys- 
tem. The twist operator is first found in the effective low 
energy model of a zigzag spin ladder 17 owing to the frus- 
trated inter-chain couplings. This parity-breaking oper- 
ator is marginal in the renormalization group (RG) sense 
and has a conformal spin l 18 . Currently, the twist opera- 
tor is believed to be the origin of incommensurate correla- 
tions in XXZ zigzag spin ladders 17,19-22 . For an isotropic 
zigzag spin ladder, the twist operator is proposed to cause 
dimerization 17 in the ferromagnetic inter-chain coupling 
region. However what role this operator plays in a spin 
tube is still unclear, which needs a further study. 



In this work, by means of the bosonization tech- 
nique and the density matrix renormalization group 
(DMRG) method we shall study a frustrated three-leg 
spin tube comprised of trigonal prism units, which derives 
from the crystal structures of [(CuC^taciiH^ClJC^ 1 and 
Ca3Co2 0g 16 , as shown in Fig. 2 (a), where the AF cou- 
plings J2 and J3 form helical paths along the tube direc- 
tion. Such a spin tube can be transformed to an equiv- 
alent spin ladder, as shown in Fig. 2 (b). Owing to the 
complex zigzag-like inter-chain couplings, the twist oper- 
ator in the bosonized Hamiltonian should depend on both 
J2 and J3, where we take the coupling J\ as an energy 
scale. Consequently, we only need to adjust J 2 and J 3 to 
study the effect of the twist term on the properties in the 
ground state. For the sake of simplicity, we consider the 
case of spin S = 1/2. Interestingly, we showed that a new 
critical phase with central charge c = § appears in this 
system, which is found in spin tubes for the first time and 
reveals a novel physical effect of the twist operator. In 
addition, one dimer phase and one conventional critical 
phase with central charge c = 1 due to the competition 
of the AF couplings J2 and J3 are also identified in the 
ground state. 



II. BOSONIZATION AND LOW-ENERGY 
EFFECTIVE THEORY 

The Hamiltonian of the present spin tube could be 
written as 

L/3 

H = ^^[Ji(S3i_2 • S 3 i_|_i + S3i_i • S3i+2 + &3i ' 83^+3) 

i=l 

+ J2{Ssi-l ■ &3i+l + &3i • S3i+i + S3; ■ S3i+ 2 ) 

+ J3(^3i-2 ■ &3i-l + S3;_i • S%i + S3i_3 ■ S3,+l)],(l) 
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FIG. 2: (Color online) (a) The structure of the frustrated 
three-leg spin tube under investigation. The circles (purple) 
indicate the sites with spin S = i>. The antiferromagnetic 
couplings J2 and J3 are presented by dashed line (blue) and 
dotted line (red), respectively. The solid line (black) indi- 
cates the AF interaction along the leg direction, (b) shows 
the equivalent spin ladder structure transformed from (a) by 
setting J3 bonds as legs. 



FIG. 1: Spin ladders and its corresponding spin tubes (the 
ladder with periodical condition along the rung direction). 

(a) shows the standard spin ladder and the typical spin tube. 

(b) and (c) are the ladders with next nearest couplings and 
the spin tubes generated by them. 
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where is the spin operator on the jth site, L is the 
total number of sites, and J\, J2, Js > arc all AF 
couplings. In what follows we consider two limit cases: 
(I) J-2 3> J3, Ji, and (II) J3 ;§> J2, Ji. In the following 
DMRG calculations J\ is set to be unity for simplicity. 



A. J 2 >J 3 ,Ji 

This case is quite simple. From Fig. 2, one may note 
that the system can be regarded as a single spin chain 
with next nearest and next next nearest neighbor inter- 
actions. The Hamiltonian can be written as 

H = J2 Si • Sj+i + J3 Si ■ Si + 2 

i i 

+ Ji^Si-Si+s, (2) 

i 

where the last two terms can be viewed as perturbations. 
Following the standard bosonization techniques 23 , in the 
continuum limit, the Hamiltonian density of Eq. (2) can 
be written as: 



g = 2 J 3 + 4 Ji - J c 



(6) 



where 3l,r are left or right moving SU(2) current op- 
erators representing the smooth magnetization part of 



spin density operators 



24 



J2<2o is the spin velocity 



with ao the lattice constant, and J c >0. Obviously, Eq. 
(3) has the same form as the Hamiltonian density for a 
zigzag spin ladder 20 except for the coefficient g of Hjj 
term. When g>0, the system is in the dimer phase with 
an energy gap: 



, const . 

A ~ exp( ). 

9 



(7) 



T-L — Hq + Hjj, 



(3) 



It is deserved to mention here that this gap could also sur- 
vive even for Ji^>J2, J3, in which the system can be seen 
as the three free spin chains with intra-chain couplings 
Ji perturbed by relevant perturbations generated by the 
inter-chain couplings J2 and J3, as in the case of spin 
ladders 24-26 . Therefore, this case as well falls into the 
gapped phase. At last, when g<0 it is in the Luttingcr 
liquid phase with central charge c—1. The phase transi- 
tion between the two phases is of a Kostcrlitz-Thoulcss 
transition 5 ' 24 . 



3 



2 




FIG. 3: (Color online) The renormalization group (RG) flows 
for the three-leg frustrated spin tube under interest with J3 
> J 2 . The variable b± = ^|±L i s the same as the definition 
of Ref. 17. The inset is the enlarged part where the RG 
flows to the strong coupling limit can be seen clearly. £ is the 
logarithmic variable in the renormalization group equations. 

B. J :i >J 2 ,Ji 

This case is very interesting. In this situation, the 
model can be reshaped into a spin ladder structure as 
shown in Fig. 2 (b). By treating all inter-chain couplings 
as perturbations, after the bosonization procedure, the 
corresponding Hamiltonian density can be written as: 

h = h + h L r + u.u + n tms t , (8) 

Vu '^y^ .: 3j, L ■ 3 jtL : + : 3 jtR ■ J,„ :), (9) 



Hlr = 9lr{3\,l ■ 3x,r + 3 2 ,l • 3 2 ,r), (10) 
Hjj = 2(Jj + J 2 )(3 1 , L + Ji,«) • (J 2 ,l + J 2 ,h), (11) 

Utwist = - 2>Jl 2 J2 ( n i^ n 2 " n 2 9 x ni), (12) 

where is the staggered part of the spin density operator 
with the subscript i the leg index, spin velocity Vs^J^a^, 
and 3j_L.R is the left or right current operator of the 
jth leg. The coefficient g^R oc — J 3 , and Hlr, Hjj and 
Htwist are marginal perturbations with scaling dimension 
d = I 17 . Htwist is the twist term that is produced by 
the frustrated inter-chain interactions. As pointed out 
in Refs. [17,19], it arouses a spin nematic phase in the 
ground state of the XXZ zigzag spin ladder. Its effect on 
an isotropic Hcisenberg model is still not very clear and 
needs further investigations. 



In fact, this present case provides a convenient way to 
study the effect of the twist term on an isotropic Heisen- 
berg system since the coefficient of Htwist could be ad- 
justed by regulating the value of J 2 . In order to look 
at the effect of a group of marginal operators, we resort 
to solve the renormalization group equations 17 ' 25 numer- 
ically. 

First, following the procedure of Rcfs. [17,24,26], 
we can rewrite Hlr, Hjj and Htwist using Majorana 
fermions (f L '^ 2 ' 3 as 

Hp er t = ai(Ci?,CiCiC£ + (l(r(r(r) 

+ a^ackaci + cicicia) 

+ ^(c^ckici + c^ciclcl 

+ clclCRCl + tlcklcl) 

+ *i(Mcici + + <&£&<2) 

+ b^cRdcUl + ckclCRCl + cUkUt), 

(13) 

where a\ = —\&2 = \a>3 = j{— 3Jl 2 ~ j2 ), 7 is a nonzero 
constant and b x = - Ji - J 2 , b 2 = + Ji + J 2 - 

The terms with the coefficient ai stem from Htwist and 
bi from both Hlr and Hjj. 

Next, we set J 2 < 3Ji to keep the coefficient of Htwist 
negative that cannot be achieved in a usual zigzag Hcisen- 
berg spin ladder. Afterwards, choosing the initial values 
of a,i and bi in the range where J3 3> J 2 , J\ , we obtain 
a new solution of the renormalization group equations 17 
as shown in Fig. 3, where the RG flows corresponding to 
Htwist go to the strong coupling regime first, which sug- 
gests a new phase may come into being. Recall that for 
a two-leg zigzag spin ladder, the RG flows corresponding 
to Hjj go to the strong coupling regime first. In the 
following we shall examine numerically the results from 
the low-energy effective theory. 

III. SPIN GAP, CORRELATIONS AND 
ENTANGLEMENT ENTROPY 

To substantiate the above bosonization analysis, we 
use the DMRG method 27 to numerically calculate the 
spin gap, spin-spin and dimer-dimer correlation func- 
tions, as well as entanglement entropy of the model. In 
the calculations, up to 8000 optimal states are kept, and 
the truncation errors are of 10~ 12 for the ground state 
and 10 -7 for the excited states. 



A. Spin gap 

Figure 4(a) is the coupling dependence of the finite-size 
spin gap at L=180 obtained in open boundary conditions 
(OBCs). With increasing J 3 for any J 2 , the spin gap has 
two area with tiny values relatively, between which there 
is a raised region of spin gap. In Fig. 4(b), J 3 depen- 
dence of spin gap at J2=2.0 is shown, where the spin 
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FIG. 4: (Color online) (a) Coupling dependence of the spin 
gap at L=180 obtained with OBCs. (b) J3 dependence of the 
spin gap for J2 = 0.5 and 2.0. Finite-size scaling of spin gap 
for (c) the gapless and (d) the gapped phases. 
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FIG. 5: (Color online) (a) The spatial dependence of the spin- 
spin correlation function (So • S3,.} in the three phases. For 
J3=0.2 and 2.0, the spins (S3 r ) are assumed along the leg di- 
rection as indicated in Fig. 2(a), while for J3=6.0, the spins 
along the leg as indicated in Fig. 2(b). (b) The spatial de- 



pendence of the dimer correlation \Dr 



phases. 
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B. Spin-spin and dimer-dimer correlation functions 

The spin-i ID gapless spin systems usually preserve 
the translation symmetry and have the spin correlations 
with power-law decay, while the gapped ones should 
break the translation symmetry and have spin correla- 
tions of exponential decay. Figure 5 (a) presents the 
spin correlation function | (So ■ S3 r )| in different phases. 
For J 2 =2.0, J 3 =0.2 and J 2 =2.0, J 3 =6.0, which are in 
the two gapless phases, | (So ■ S3 r )| decays with a power 
law, while at J2=J 3 =2.0 that is in the gapped phase, 
| (So ■ S 3r )| has an exponential decay with a short corre- 
lation length £ ~ 2.7. The breaking of translational sym- 
metry is usually detected by the dimer-dimer correlation 
functionI)( i) j) ) ( fcj j)=((Si-S J ')(Sife < Si))-(Sf-S i7 ')(Sife'Si). In 
Fig. 5(b), |-D(o,3),(3r,3r+3)| is shown for the three phases. 
In both gapless phases, the dimer correlations decay with 
a power law, which is consistent with the preserved trans- 
lation symmetry and gapless spin excitations. In the 
gapped phase, the dimer correlation builds a long-range 
order (LRO). 



C. Entanglement entropy 



gaps increase to the maximum value in the raised region 
and decrease slowly with growing J 3 . while those in other 
regions diminish fast, probably suggesting a gapped spin 
excitation in the middle region and a gapless excitation in 
other regimes. As shown in Figs. 4(c) and (d), the spin 
gaps away from the raised region extrapolate to zero, 
while those within the raised region go to finite values in 
the thermodynamic limit, therefore indicating two gap- 
less and one gapped phases. These results are also in 
accordance with our previous bosonization analysis that 
the system will be in the gapless phase for J2^J3, J\ and 
J2, Ji«CJ3, and is gapped in other parameter regimes. 



To distinguish the two gapless phases 28 , we calculated 
the central charge c from the block entanglement entropy 
of the system Si=— Tr/?/ lnp;, where pi is the reduced 
density matrix of a subsystem of size /. For a gapless 
system with periodic boundary conditions (PBCs), the 
entanglement entropy is given by 



Si = - In 



~N . , 




— sin 




7T ' 





ffPBCs , 



(14) 



where N is the total length of the system, and gpBCs is a 
nonuniversal constant with PBCs 29 . The central charge 
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FIG. 6: (Color online) Fitting of the entanglement entropy 
given by Eq. (14) to the DMRG results in the two gapless 
phases. The central charges from the fitting are obtained to 
be c=l and 3/2, respectively. 



is obtained by fitting the entanglement entropy given by 
Eq. (14) to the DMRG results for the system with L up 
to 360. 

As shown in Fig. 6, in the gapless phase with large 
J 2 (for instance J 2 =2.0, J 3 =0.1), both results for Si are 
fitted quite well with c=l, which indicates that this gap- 
less phase belongs to the same universality class as the 
spin-i Heisenberg AF chain. In the gapless phase with 
large J3 like J2=0.2, Js=7.0, the central charge is identi- 
fied as c=|. Therefore, we can label the different phases 
with its central charge in Fig. 4 (a), which comprises the 
phase diagram. As the model preserves the SU(2) sym- 
metry, the transition from the dimer phase to the gapless 
phase with c=| might be in the SU(2)fc = 2 Wess-Zumino- 
Witten-Novikov (WZWN) universality class 23 . In com- 
bination with the bosonization analysis, it is observed 
that for Jz^>J2 the RG flows of Htwist g° to the strong 
coupling regime faster than Hjj that leads to a dimer- 



ized phase 17,24 , and one may judge that such a gapless 
phase with nontrivial central charge c=| in the frustrated 
three-leg Heisenberg spin tube is probably induced by the 
twist term Htwist- 



IV. SUMMARY 

To summarize, by means of the bosonization technique 
we develop a low-energy effective theory for the spin- 1/2 
frustrated three-leg spin tube in two limit cases, and also 
invoke the DMRG calculations on the spin gap, spin-spin 
and dimcr-dimer correlation functions as well as the en- 
tanglement entropy to elucidate the effective analyses. 
We have discovered a dimer phase and two distinct criti- 
cal phases in this system. The dimer phase is character- 
ized by the finite spin gap, exponentially decaying spin- 
spin correlations, and a dimer-dimer LRO. The critical 
phases are found to have gapless spin excitations, power- 
law decaying spin and dimer correlations. The differ- 
ent central charges 1 and I distinguish the two critical 
phases. Based on the bosonization analysis, the novel 
critical phase with c=| can be attributed to the nega- 
tive sign of Htwist in Eq. (12), which in turn reflects a 
new effect of the twist term, and might belong to the 
SU(2) fc=2 WZWN universality class. The RG flows in 
this frustrated three-leg Heisenberg spin tube with J 3 3> 
J 2 differ from those of the two-leg zigzag spin ladder. 
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